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Abstract 

We calculate the specific heat of the antiferromagnetic spin-1 chain com- 
pound Y2BaNii_a;Zn2.05 in the presence of a magnetic field. The low-energy 
spectrum of a Heisenberg Hamiltonian, which includes realistic anisotropics, 
has been solved using the density matrix renormalization group. The ob- 
served Schottky anomaly is very well described by this theory. For large 
chains {N > 50), contrary to previous interpretations of the specific heat 
data, we find S = 1/2 states. These results are thus also consistent with 

electron-spin-resonance data for NENP. 
PACS numbers: 75.40 Cx, 75.10Jm, 75.40 Mg. 
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A great deal of interest in one dimensional Heisenberg chains with nearest-neighbor an- 
tiferromagnetic (AF) exchange coupling, J, has been originated by Haldane's conjecture 
that integer-valued spin chains would exhibit a gap in the spin-wave excitation spectrum 

in contrast to half- integer spin chains which would be gapless with a linear dispersion 
relation above the ground state. This quantum many-body phenomena is different from the 
usual source of gaps in magnets, namely, single-ion anisotropy, which does not involve cor- 
relation effects. Subsequent observations of large gaps in spin-1 quasi-one-dimensional sys- 
tems present in the structure of some compounds like CsNiCls 0, [Ni(C2HgN2)2(N02)]C104 
(NENP) [^^, and Y2BaNi05 [§,0, have confirmed Haldane's observation. 

Affleck et al. have shown that the exact ground state of the Hamiltonian Sj • 
Sj+i + (Sj ■ Sj+i)^/3 is a valence-bond-solid (VBS) state. In this state, each S =1 spin 
is represented by symmetrization of two S = 1/2 entities. These S = 1/2 spins at each 
site are coupled with nearest-neighbor S = 1/2 spins, one to the left and the other to the 
right, to form singlets. Hence, an open chain has two unpaired 5=1/2 spins, one at each 
end. They have also proposed that the VBS state describes qualitatively the physics 
of the Heisenberg model (without the biquadratic term). In agreement with this proposal, 
exact diagonalization |Q| of finite open chains shows that the four lowest-lying states are a 
triplet and singlet whose energy separation approaches zero exponentially with increasing 



length. Monte Carlo and density-matrix renormalization-group (DMRG) |TT[| studies 



clearly show the presence of S = 1/2 end states. This picture is also supported by electron 



paramagnetic resonance (EPR) measurements of NENP doped with non-magnetic ions ||T2 
where resonances corresponding to the fractional spin S = 1/2 states at the "open" ends of 
the Ni chains were observed. 

However, Ramirez et al. also tested the presence of free 5 = 1/2 states by studying 
the specific heat of defects in NENP and Y2BaNi05, with magnetic fields up to 9T and 
temperatures down to 0.2K. On the basis of a comparison between two very simplified 
pictures, they found that the shape and magnitude of the Schottky anomaly associated 
with the defects in Y2BaNii_2:Zna;05 are well described by a simple model involving spin-1 



excitations, instead of the S = 1/2 excitations of the VBS. These results are in contrast 



with the above mentioned EPR measurements of NENP |jT2| . 

While previous theoretical studies of the spin-1 Heisenberg chain confirm the 

existence of free 5 = 1/2 spins at the ends of sufficiently long chains in agreement with 
the EPR measurements, it seems difficult to reconcile these results with the specific heat 
measurements and no theoretical calculation of the latter exists so far. 

In this Letter, we solve the low energy spectrum of the appropriate Hamiltonian which 
describes the Ni chains of Y2BaNi05 using DMRG and finite size scaling. Assuming that 
defects are randomly distributed (which leads to a Poisson distribution for the length of the 
open chains), we can calculate the specific heat for any concentration of non-magnetic im- 



purities. Our results show that there is no contradiction between EPR |12] and specific heat 



measurements in Y2BaNii„^Zna,05 |T^. They confirm the existence of spin-| excitations for 



sufficiently long chains and are in excellent agreement with the curves measured by Ramirez 
et al. ||13[. 

Y2BaNi05 has an orthorhombic crystal structure with the Ni^+ {S =1) ions arranged 
in linear chains with a nearest-neighbor AF superexchange coupling. The interchain cou- 
pling is more than three orders of magnitude weaker, making this compound an ideal one- 
dimensional antiferromagnetic chain. While each Ni atom is surrounded by six O atoms in 
near octahedron coordination, the true site symmetry is D2hi and the appropriate Hamilto- 



nian reads 14-16 



i 

where z is along the chain axis and Stot is the total spin. Recent estimates based on fits 
of the Haldane gaps (in x, y, and z directions) measured by inelastic neutron scattering , 

indicate J ~ 280/^ , D ~ -.038 J , and E .013 J |4|,|15l. 

As pure Y2BaNi05 has a Haldane gap ~ lOOi^ []7| JT^ , p!5[] , the spin wave contribution to the 



specific heat C{T) is negligible below 7K. In this temperature range, C(T) is dominated by 
the effect of the defects, which is manifested as a 1/T^ rise in C(T), and a Schottky anomaly 



in the presence of an applied magnetic field. The height and width of the Schottky anomaly 
strongly depends on the spin value of the low energy excitations. For these reasons, it is 
necessary to solve the low energy spectra {uj < lOK) of H for all values of to have an 
accurate theoretical picture. 

By means of the DMRG method, we have calculated the two lower eigenenergies in 
the S^^^ = subspace and the lower one with S^^^ = 1 (this state is degenerated with the 
^tot — ~^ due to the time-reversal symmetry of H) for all < 40 and E = 0. The energy 
difference between the excited states and the ground state decays exponentially to zero with 
increasing A^. This behavior, shown in Fig. 1 for an arbitrary value of anisotropy D = —0.1, 
allows us to extrapolate the energies to all values of > 40, and demonstrate that the 
free S =1/2 spins at the end of the chain persist in presence of anisotropy. This issue is 
easy to understand considering that S =1/2 spins cannot be affected by anisotropy due to 
time-reversal symmetry. At this level, it is important to remark that the two S =1/2 spin 
excitations, predicted by VBS, have a finite localization length / ~ 6 sites p|-pJ]|. Therefore, 



while they are nearly free for large (A^ >> I) open chains, the interaction between them is 
considerable when the length of the chain A^ is comparable to 21. This interaction splits the 
two S =1/2 states into a singlet (ground state for even A^) and a triplet (ground state for 



odd A^) one [|T7 



The difference between any two energies of the above mentioned low-energy states is 
linear in D, and the quadratic corrections are negligible |T^. This widely justifies the 



validity of perturbation theory to first order in D. By symmetry then we can include the 
term J^iEKSf)"^ — (Sf)"^] to first order. Thus we find the following low-energy effective 
Hamiltonian including a spin S = 1 and a singlet state |0): 

Heff = (EoiN) + a{N)D) |0) (0| + Df3{N)S', 
+E(3{N){Sl-S'y)-gf,BB-S 

where Eq{N), a{N), and f3{N) are functions of the chain length A^ (determined from the 
DMRG data). The validity of the last term has been verified explicity by calculating the 



matrix elements of 5'^^, and 5'^^ for all chains. Heff determines the thermodynamics of the 
system at temperatures well below the Haldane gap. 

For a random distribution of defects, the specific heat per elementary cell of Zn-doped 
YsBaNiOs is: 

oo 
N=l 

where Cn{B,T) is the specific heat of a segment of length described by -ffe//, and x is 
the concentration of missing Ni atoms. 

In order to minimize the effect of the lattice on C(T), Ramirez et al. ||13| plot the data, 
measured in Y2BaNio.96Zno.04O5 , as the difference [C(3T) — C(6T)]/T . In Fig. 2, we 
compare this difference calculated with -ffe//, with the data taken from Ref. |]13|. As the 



samples were powdered in order to ensure good thermal equilibration [0], the theoretical 



curve was obtained by averaging over all magnetic field directions. We have taken J = 280K 



, D = — 0.038J consistent with neutron scattering experiments |T^JT5[] and E = — 0.032J, 
somewhat larger than the experimental value. We have also considered g = 2.35 which is 
a reasonable value for the g factor [|T2],|TB|. The agreement between experiment and theory 
is excellent. Note that although the fit with the singlet-triplet model is good, the result we 
obtain here is better, particulary for temperatures around 0.5K, and between 3 and 6K. We 
also used a doping concentration of x = 0.0428 higher than the nominal value x = 0.04. This 
is due to the existence of native defects, also present in the pure compound. By fitting the 
data corresponding to pure Y2BaNi05 we obtain x = 0.0028 instead of the value x = 0.008 
obtained with the singlet-triplet model in Ref. [jl3[ (see Fig. 3). Note that for small values 
of X, the height of the peak is proportional to x. 

From Figs. 2 and 3 it is clear that the inclusion of interactions between end S =1/2 
spins through the chain via the full Hamiltonian H, together with anisotropy, eliminates the 
apparent discrepancy between VBS theory predictions and measured specific heat [|13|. In 
Fig. 4 we have separated the contributions of different chain lengths to [C(3T) — C(6T)]/T . 
The full line corresponds to the result for non-interacting S =1/2 end spins. This behavior 



is obtained for ^ 50 where the singlet-triplet gap is much smaller than q^bB. The 
effect of anisotropy is stronger for shorter chains, changing the position and height of the 
Schottky peak. This behavior is mainly due to odd chains because they have an S =1 
ground state. The anisotropy splits the 5*2 = and Sz = ±1 (ground state for D < ) 
triplet states, and the energy separation A(A^) decreases with (as can be inferred from 
Fig. 1). As one can see from Fig. 4, for large chains (A^ > 30), the result approaches 
asymptotically the VBS prediction, i.e. the peak is higher and shifted to slightly smaller 
temperatures as compared to the experimental result. However, the curves for shorter chains 
(A(A^) ^ q^bB) show the opposite behavior, and both tendencies are compensated in the 
final result. To understand this effect let's consider the simplest case of E = = By = 
and odd A^ . The energy difference corresponding to the Schottky anomaly, E{Sz = 0) — 
E{Sz = —1), becomes q^bBz + A(A^) when a magnetic field is applied. This explains the 
shift to the right registered for the shorter chains. As a consequence, the sum over small 
(A^ ^ 20) values of A^, weighted with the Poisson distribution, yields a Schottky peak of 
larger width and lower height . The anomalous features at low temperature in Figs. 2 and 
4 are due to the contributions of chains with even A^ for which the difference between the 
energy of the triplets with = ±1 and the singlet ground state becomes of the order of 
the Zeeman energy g^sB. For 6T this crossing occurs for A^ between 18 and 20 giving rise 
to the large negative contribution for the corresponding curve in Fig. 4. For B = 3T the 
crossing occurs for A^ between 22 and 24 (see the 21 < A^ < 30 curve in the same figure). 

Our results allow us to understand why the S = 1 states of short chains were not seen 
g in the EPR experiments in NENP fT^. For D/J ~ 0.2 the splitting A(A^) is greater 



or of the order of the applied magnetic field gfiBBa {Ba ~ 0.35T ) when A^ < 38, while for 
A^ > 54, A(A^) < O.lgfisBa ■ These larger chains (76% of the total number of chains for 
X = 0.005 ) contribute to the same EPR signal with total weight 0.0038 per transition metal 
atom. On the other hand, each of the chains with A^ < 38 ( weight ^ .0002), give place to 
a different gyromagnetic frequency larger than the one they considered . 



In conclusion, by solving the low-energy spectrum of a Heisenberg Hamiltonian H which 



includes axial and planar anisotropy we have reproduced the low-temperature and high- 
field specific-heat data measured in Y2BaNii_^.Zna,05 |T^. Our results are consistent with 
valence-bond-solid predictions of S =1/2 end chain excitations, which are asymptotically 
free for large chain segments. However, their interaction for short segments are critical for the 
understanding of the specific heat behavior. These results remove the apparent discrepancy 
between the specific heat data for Y2BaNi05, interpreted in Ref. |13| with a singlet-triplet 
model, and electron paramagnetic resonance data for NENP. 
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Figure 1: Difference between the energy of the triplet = ±1 {Sz = 0) and the singlet 
states in the presence of anisotropy D are denoted with circles (squares) . a) even, b) 
odd. The full line is an exponential fit of these differences for > 19. 



Figure 2: Specific heat difference between 3 and 6T for Y2BaNio.96Zno.04O5 []T3|. The 
dashed and solid lines correspond to fits with the prediction oi H (J = 280K, D = 
-0.038 J, E = -.032 J, X = 0.0428, g = 2.35) and the singlet-triplet model (x = 0.048). 

Figure 3: Specific heat difference between 3 and 6T for undoped Y2BaNi05. Solid circles 



are experimental data ||13[. The solid line corresponds to a fit with the prediction of H 
for the same parameters as in Fig. 2 except for x = 0.0028, and the dashed line to the 
singlet-triplet model with x = 0.008. 

Figure 4: Contributions of chains in different length intervals to the specific heat differ- 
ence between 3 and 6T calculated with H for the same parameters as Fig. 2. The full line 
corresponds to the prediction of VBS theory for non-interacting S* =1/2 moments. The final 
result is obtained adding each curve with the weight denoted inside the figure. 
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Fig. 1 Hallberg et. al. 
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Fig. 3 Hallberg et. al. 
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Fig. 4 Hallberg et. al. 



